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ABSTRACT 

Magnetars are known to have optical and/or infrared emission, but the origin of the emission is not 
well understood. In order to fully study their emission properties, we have carried out for the first time 
optical linear polarimetry of the magnetar 4U 0142+61, which has been determined from different 
observations to have a complicated broad-band spectrum over optical and infrared wavelengths. Erom 
our /-band imaging polarimetric observation, conducted with the 8.2-m Subaru telescope, we deter¬ 
mine the degree of linear polarization P = 1.0+3.4%, or P <5.6% (90% confidence level). Considering 
models suggested for optical emission from magnetars, we discuss the implications of our result. The 
upper limit measurement indicates that different from radio pulsars, magnetars probably would not 
have strongly polarized optical emission if the emission arises from their magnetosphere as suggested. 

Subject headings: X-rays: stars — stars: neutron — pulsars: individual (4U0142+61) — polarization 


1. INTRODUCTION 

It is currently accepted that magnetars are highly mag¬ 
netized neutron stars whose dipole magnetic field reaches 
~10^^-10^® G, implied by their slow spin periods of 2- 
12 s and rapid spin-down ra tes of 10~^^ — 10~^° s s~^ 


(see reviews, e.g., given bv IWoods & Thomuson 

2006 

Kasoil 120071 IMereghettil 120131 but see iRea et al. 

2010 


20121 for two low magnetic field magnetar cases recently 


discovered). Neutron stars which have traditionally been 
classified as anomalous X-ray pulsars (AXPs) or soft 
Gamma-ray repeaters (SGRs) belong to this magnetar 
class, and over 20 magnetars have been identified up to 
now (jOlausen fc Kasnl 120141) . Their bright X-ray lumi¬ 
nosities cannot be explained by rotational energy loss 
rates and hence the energy source is considered to be 
provided by the decay of their ultra-strong magnetic 
fields. AXPs and SGRs are well k nown to show a vari¬ 
ety o f high-energy phenomena fe.g.. lWoods fc ThompsonI 
120061 : lKasDil[200% including the exceptionally bright ‘gi¬ 
ant flares’ wh ose peak luminosity amounts to 10"^^ 
erg s~^ fe.g.. iHurlev et aH 119991 iTerasawa et al.l 120051 : 
iTanaka et al.ll200^ ^ 

In addition to these intensive studies at X-ray ener¬ 
gies, magnetars have been subjected to multi-wavelength 
follow-up observations in radio, optical, and infrared 
(IR) bands. Among the magnetars identified so far, 
the AXP 4U 0142+61 is the best studied magnetar at 
optical and IR wavelengths due to its relatively short 
distance (di stance d ~ 3.6 kpc) and low extinction 
(Ay — 3.5; iDurant fc van KerkwiikI l2006allbll An op¬ 
tical c ounterpart was first discovered bv IHulleman et ^ 
(j2000fl from 4U 0142+61. The optical flux was higher 
than the Rayleigh-Jeans tail of the blackbody component 
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seen in X-ray band, requiring another emis sion niecha- 
nism. Subsequently, IKern fc Martini (1200211 found that 
the optical emission is pulsed at its spin period and 
that the pulsed fraction is ~27%, whic h is higher than 
4%-1 4% observed in its X-ray emission (|Gonzalez et al.l 
1201011 . A very faint near - IR co unterpart was also iden¬ 
tified by iHulleman et al.l (1200411 , and then Spitzer Space 
Telescope detected mid-IR emission from 4U 0142+61 
(|Wang. Ghakrabart v . fc K aplanll2006ll . 

Importantly. iW^g et aLh 20061) found that two differ¬ 
ent emission components are needed to explain the com¬ 
bined optical and IR spectral energy distribution (SED). 
Namely, the thermal blackbody-like component, possibly 
emitted from a debris disk around the pulsar, is domi¬ 
nant over the 2.2-8 /im range, while another power-law 
like component (represented by oc can account 
for the re sidual emission in t he optical VRI and near-IR 
J bands (|Wang et al.ll2006ll . Follow-up Spitzer mid-IR 
spectroscopy and 24 fim imaging of the AXP were con¬ 
sistent with the interpretation of a debris disk for the 
IR emission (|Wang et al.ll200§l . However, the origin of 
the optical power-law emission is still unclear, although 
it would probably arise from the magnetosphere because 
the emission is relatively highly pulsed. 

In order to shed new light on the optical emission 
mechanism of magnetars, we have carried out polari¬ 
metric observations for the nearest and least absorbed 
m agnetar 4U 0142+61 using the Su baru 8.2-m telescope. 
In [Wang. Tanaka, fc ZhongJ (|2012f) . we have reported a 
4.3% upper limit (90% confidence level) on the degree of 
circular polarization in /-band emission from the source. 
In this paper we present the result from our linear po¬ 
larimetry of the source at /-band. 

2. OBSERVATION AND DATA REDUCTION 


We carried out linear imaging polarimetry of 4U 
0142+61 at I-band with the 8.2-m Subaru Telescope on 
2013 December 22-23. Two first half of the telescope 
time in each night were awarded to us by the National 
Astronomical Observatory of Japan. Due to high humid¬ 
ity in the first night, only three hours of data were taken. 
For the second night, full five hours of the telescope time 
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were used. 

The Faint Object Ca mera and Spectrograph (FOCAS; 
iKashikawa et al.ll2002f) was used for the polarimetric ob¬ 
servation. In this mode, a Wollaston prism and a half¬ 
wave retarder are inserted to the collimated beam. An 
incident beam is splited by the Wollaston prism into two 
orthogonally polarized beams, one ordinary (o-beam) 
and the other extraordinary (e-beam). To avoid blend¬ 
ing of the two beams, a standard focal mask has to be 
used. In order to measure the degree of linear polariza¬ 
tion and the position angle of polarization in the celestial 
plane, one set of four exposures of a target with the half¬ 
wave plate at four position angles, 0?0, 45?0, 22?5, and 
67?5, are required. The FOCAS detector was two fully- 
depleted-type 2kx4k CCDs, the pixel scale of which was 
0'.'104pixel“^. We 2x2 binned the detector for our ob¬ 
servation. 

In our observation, each set of four exposures were 
taken by setting the half-wave plate at the four posi¬ 
tion angles alternately. In order to avoid possible sever 
saturation caused by bright stars, the time of each expo¬ 
sure was 2 min. Between the sets of the exposures, the 
telescope was five-point dithered to avoid bad pixels on 
the CCDs and to help remove cosmic ray hits during the 
data reduction. In total, in the first and second night 
we obtained 19 sets and 32 sets of the exposures, re¬ 
spectively. The observing conditions were unfortunately 
mediocre comparing to what we wished (075 seeing would 
be ideal). The average seeing was approximately 0780 in 
the first night and 0786 in the second night. 


• • 


Upper panel 


• • 

Bottom panel 


Fig. 1. — Subaru/FOCAS linear polarization image of the 
4U 0142-1-61 field at I-band. The two polarized beams were 
recorded at the upper and bottom panels, in which the AXP is 
indicated by X. Several ghost stars, marked by dotted circles, are 
present in the bottom panel. 

In order to determine the zero point of the position 
angle of polarization in the celestial plane, the standard 
star BD-|-64dl06, which has strongly polarized emission, 
was observed. 

We used the IRAF packages for data reduction. The 
images were bias subtracted and flat fielded. Dome flats 
at the four position angles of the half-wave plate were 
taken and used for flat fielding respectively. For our tar¬ 
get, the images made at each angle were then combined 
into one final image of the target field by positionally 
calibrating them to a reference image. In the first night 
at the position angle of 45?0, the target position in one 
image was contaminated by a cosmic ray hit, and the 
image was not included in image combining. In total 



the on-source exposure times at each position angle is 38 
min (36 min for the position angle of 45?0) and 64 min in 
the first and second night, respectively. Because of the 
faintness of the target and mediocre seeing during our 
observations, the signal-to-noise ratios of the target in 
each night’s images are not sufficiently high. We there¬ 
fore combined all the data of the two nights and made 
a set of four images of the target field at the four posi¬ 
tion angles respectively. A target-field image, made by 
combing all the images, is shown in Figure [TJ 
For the standard star, aperture photometry was per¬ 
formed and the o-beam and e-beam brightnesses at each 
position angle were obtained. For the target, because of 
its faintness, we used the point spread function (PSF) 
photometry tasks from IRAF’s DAOPHOT package to 
measure its brightnesses. We also derived aperture cor¬ 
rections by using 8 in-field, relatively bright stars (bright¬ 
est stars in the field were saturated). A radius of 10.0 
pixels (~ 271) was used for aperture photometry of these 
bright stars. The uncertainties on the corrections are 
negligible, as they are much smaller than those from pho¬ 
tometry of the target. These magnitude measurements 
and uncertainties are given in Table [TJ 


3. RESULTS 

The degree of linear polarization P is calculated from 


P=^Q^ + U^ , ( 1 ) 


where Q and U are Stokes parameters. They are deter¬ 
mined from 


Q = 


1 — ai 
1 + fli 


U = 


1—02 
1 -b 02 ’ 


where 


„ _ , '* 0?0 „ _ r^22°5 

Ol — W-, 02— - 


Here k = le/Io, the ratio of the intensity I of the target 
in the e-beam and o-beam frames at each of the four 
position angles of the half-wave plate. The position angle 
ip of linear polarization is determined from 


^ = iarctan(^) . 


( 2 ) 


Then from standard propagation of uncertainty, where a 
first-order Taylor series expansion is used and the vari¬ 
ables in the functions are assumed to be independent, the 
uncertainties on P and tp can be calculated from those o n 
intensity measurments (see, e.g.. iKawabata et al.lll999ll . 

Using the formulae given above, we first checked the 
polarization measurements for the standard star. From 
the intensity measurements, we found P = 4.710 ± 
0.046% (Q = 2.325 ± 0.046%, U = -4.096 ± 0.046%) 
at /-band, whi ch is consistent with the report ed value 
4.696±0.052% (iSchmidt. Elston. &: LuDidll992h . There¬ 
fore no instrumental correction to the polarization mea¬ 
surement was needed. The position angle tp' in the in¬ 
strumental coordinate (cf. Equation (2)) was found to 
be tp' = —30?21 ± 0?28. Comparing it to the reported 
value tp = 96?89 ± 0?32 (iSchrnidt et al.lll9^ . which is 
the position angle in the celetial plane, the correction S 
was determined to he S = tp' — tp = —127710 ± 0743. 
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We then calculated the Stokes parameters and degree 
of linear polarization for 4U 0142+61, and found Q = 
-0.97 ± 3.40%, U = -0.21 ± 3.31%, and P = 1.0 ± 
3.4%. Since the uncertainty dominates, a 90% confidence 
level constraint on the degree of linear polarization was 
derived. The upper limit is 5.6%. In addition, also due 
to the large uncertainties on U and Q, ip = tp' — 5 = 
43° ± 96° (where ip' = —83?9), which provides no useful 
information. 


4. DISCUSSION 

Neutron stars are known to generally have faint op¬ 
tical emission, either the Rayleigh-Jeans tail of high- 
temperature thermal surface emission or nonthermal 
emission arising from their magnetospheres. Optical 
(or near-lR) fluxes of magnetars lie above the spec¬ 
tra extrapolat ed from their blackbody-like X-ra,Y com - 
ponents (e.g., iHulleman et al.l [20001 : iMereghettil 1201311 . 
which excludes the former case for magnetars. For the 
latter, polarized emission is expected and actually is 
observed in optical emi ssion from radio pulsars (e.g., 
iSlowikowska et all 12009( 1. While the current emission 
models, such as the polar cap, outer gap, two-pole caus¬ 
tic, an d striped pulsar wind ( see the detailed discus¬ 
sion in ISlowikowska et all 120091 and references therein), 
which involve radiation mechanisms including curvature, 
synchrotron, or inverse Compton scattering radiation, 
can not fully explain the well-studied Crab pulsar case, 
5-10% phase-averaged linear polarization is detected 
in optical emission from seve ral close or young radio 
pulsars (e.g., the Crab pulsar. 


Moran et ahf [Ml B0 540-69 


Slowikowska et al.l 


Middleditch et al.l 


200 


198 


Lundnvist et al.l 20111: t he Vela puls ar, Mignani et ahl 
20071: iMoran et all 12014 Rn65 6-H4. IKern et al l 12001 
Mignani et al.ll20ll B1509— 58. lWagner fc Seifertll2n0fl[l . 


The polarization upper limit of 5.6% we have obtained 
suggests, while marginally, a different emission mecha¬ 
nism from that considered in radio pulsars and thus sup¬ 
ports a class of their own for magnetars. 

Considering the surface magnetic fields of R ~ 10^® G 
for magnetars, lEichler et al.l (j2002l l have suggested that 
optical emission from them could be similar to radio 
emission from radio pulsars, i.e., due to synchrotron ra¬ 
diation from electron/positron pairs. In this scenario, 
strong linear polarization should be seen, since pul¬ 
sars’ radio emission is known to have the degrees of 
linear polarization (phase-average d) in a range of from 
10 % to as high as 100 % (e . g.. ICould fc Lvnel 119981 : 
IWeltevrede fc JohnstonI |200^ iHan et al. 200911 . Ad- 
ditionallv. [Beloborodov fc Thompson! ( 200711 have sug¬ 
gested ion cyclotron emission or curvature emission by 
electron/positron pairs as two possible mechanisms for 
magnetars’ optical emission. For the both mechanisms, 
certain degree of linear polarization might be expected. 
Our measurement suggests zero or low linear polarization 
in optical emission, not supporting the scenarios. 

However, the propagation effects in a magnetosphere 
may cause a strong depolarization of optical emission. 
The natural wave modes of optical wave in the pulsar 
magnetosphere are usually two orthogonal linearly po¬ 
larized modes: the ordinary mode polarized in fc — R 
plane and the extraordinary mode perpendicular to that 
plane. The adiabatic evolution condition of the two lin¬ 


ear modes is (see details in iWang. Lai, fc Hanll2010l l 

~ 4 X 10 - 577373 - 81 *^ 1 - 5 '* » 1- (3) 

Here rj = N/Nqj is the multiplicity, = 77 / 10 ^, 73 = 
7/195 with 7 the Lorentz factor of the streaming plasma, 
R 5 = R/IO^G and 1^15 = iz/lO^^Hz are the magnetic field 
strength and wave frequency. Note that we reasonably 
suppose (p'g ^ 1/ric (ric is the light cylinder radius). As¬ 
suming a dipole, the magnetic field strength can be writ¬ 
ten as 

R = R,i5(r/R,)-5 G = 9.2R,i5R-5(r/rie)-5 G, (4) 

where R* is the surface magnetic field, R*i 5 = R/lQi^G, 
R* the neutron star radius, and we set R* = 10 km, 

Pios = R/lOs. For r ^ we 

have I Fad I "C 1. Thus if it is in the outer magnetosphere 
where magnetar optical emission is generated, the mode 
evolution would be usually non-adiabatic, implying that 
the polarization direction almost does not change when 
propagating through the magnetosphere. 

Gonsidering the emission of a relativistic particle ac¬ 
celerated along the magnetic field direction, the initial 
polarization state should be the pure ordinary mode 
{E II Rj_) in each direction of the I /7 emission cone. 
The polarization percentage of the combined wave of 
the cone obvi ously equals to zero a t the e mission point. 
As pointed bv IGheng fc RudermanI (Il979f) . the polariza¬ 
tion directions of the cone could be aligned by adiabatic 
walking after propagating a distance, which gives a very 
high polarization percentage such as in radio emission. 
However, if the optical emission comes from the outer 
magnetosphere, adiabatic walking would not occur since 
the mode evolution is non-adiabatic. The final polariza¬ 
tion degree should be close to zero in this case. In this 
emission model, if the emission region extends to inner 
magnetosphere where mode evolution is adiabatic, par¬ 
tial polarization are expected, which may be the case 
of the Crab pulsar. For 4U 0142+61, the optical emis¬ 
sion region could be higher so that the mode evolution 
is purely non-adiabatic everywhere, and the final polar¬ 
ization degree is close to zero. 

Even the initial emission at each height is highly po¬ 
larized for some other mechanisms, the total emission 
from the whole magnetosphere could also be depolar¬ 
ized due to the aberration/retardation effect. The ob¬ 
served emission at one phase may come from different 
magnetic field curve planes where the local polarization 
directions are different. Different from the radio band, 
pulsar optical emission may extend from the inner mag¬ 
netosphere to light cylinder radius, such as the whole 
slot gap and outer gap region, which makes the depo¬ 
larization due to the aberration/retardation effect much 
stronger. The detailed depolarization degree depends on 
the emission geometry in a magnetosphere. We also note 
that possible precession of 4U 0142+61 with a pe riod 
of ~15 hours was reported (jMakishima et all 120141 1. al- 
though it was not con firmed by NuSTAR observations 
([Tendulkar et al.ir2015ll . This precession could cause the 
polarization degree to be further averaged down in our 3 
and 5 hours observations. 

In addition to the pulsar magnetospheric origin con- 
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sidered above, optical I-band emission from 4U 0142+61 
could have other origin. For example, part of it (the 
non-pulsed component) could arise from the debris disk 
around the ma gnetar. In a recently re-proposed model 
for magnetars (iMalheiro et al.l l2012t following the pi(3 - 
neering work bv iMorini et al.lll988l an d iPaczvnski Ill99flll . 
they are suggested to be massive, highly magnetized 
white dwarfs with emission powered by their rotational 
energy (similar to pulsars). The optical and IR SED 
of 4U 0142+61 then might be explained by the thermal 
surface emission from a massive white dw arf plus that 
from a surrounding disk (for details, see iRueda et al.l 
l2013f) . For these possibilities, low linear polarization (at 
most a few percent) would be expected, since the central 
star’s emission can significantly lower the polarization 
level in optical light from a debris disk system (although 
emis sion from a debris disk may have high polarization; 
e.g.. lGarcia fc G6mezll201^ . or thermal sur face emission 
from a star cannot be highly polarized ('e.g.. lCheng et al.l 
I1988I1 . Our constraint on the linear polarization is actu¬ 
ally consistent with these possibilities. 

Our Subaru polarimetry likely represents the best ef¬ 
fort for measuring polarization in optical light from 
4U 0142+61 or the known magnetars in general, since 
they do not ha ve detectable optical e mission or are 
extremely faint (lOlausen fc Kaspil 1201411 . However, it 
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would be interesting to carry out such observations when 
they, particularly 4U 0142+61, are in an outburst. The 
comparison would provide further information for our 
understanding of magnetars’ emission mechanisms and 
related properties. In the future, polarimetry with an 
extremely large telescope, such as the Thirty-Meter Tele¬ 
scope, would certainly help our understanding by ob¬ 
taining a much tight constraint or a measurement in a 
much shorter observation. The result would possibly help 
determine the emission mechanism at optical bands for 
magnetars and even identify the true nature of magne¬ 
tars. 
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TABLE 1 

PSF FITTING PHOTOMETRY OF 4U 0142+61 


Sub-image 

^psf 

(mag) 

(mag) 

fTlr — lO 
(mag) 

e-beamQOQ 

26.809+0.071 

0.250 

26.559+0.071 

o-beamQOQ 

26.858+0.072 

0.246 

26.612+0.072 

e-beam^gOp 

26.780+0.073 

0.255 

26.525+0.073 

o-beam^gOp 

26.792+0.078 

0.256 

26.536+0.078 

e-beam220 5 

26.693+0.065 

0.257 

26.436+0.065 

o-beam220 5 

26.846+0.081 

0.253 

26.593+0.081 

e-beamg,j,Og 

26.731+0.060 

0.260 

26.471+0.060 

o-bcarn^-o p. 
_ 67 . 5 

26.874+0.076 

0.255 

26.619+0.076 


Note: instrumental magnitude m = 28 — 2.5log(fiux); rrip^f, Arricor, and are magnitudes obtained from PSF fitting, aperture 

corrections to a radius of 10 pixels (the uncertainties are negligible), and corrected magnitudes. 



